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Abstract 

Effects from a finite top quark mass on differential distributions in the Higgs+jet 
production cross section through gluon fusion are studied at next-to-leading order 
in the strong coupling, i.e. 0(a^). Terms formally subleading in l/m t are calculated, 
and their influence on the transverse momentum and rapidity distribution of the Higgs 
boson are evaluated. We find that, for the differential K-factor, the heavy-top limit 
is valid at the 2-3% level as long as the transverse momentum of the Higgs remains 
below about 150 GeV. 

1 Introduction 

The past few years in particle physics have been characterized by ever more sensitive 
exclusion limits for Higgs bosons (see, e.g., Refs. [l]-[3] ) . These results were based on the 
combination of experimental data, background calculations and extrapolations, theoretical 
expectations for the signal, and careful estimates of the associated uncertainties [4j[5]. 
Concerning the signal cross sections, these uncertainties have several sources: the parton 
densities (PDFs), the strong coupling a s (Mz), and higher order perturbative effects, for 
example. 



An uncertainty which is very specific to hadronic Higgs production in the Standard Model 
concerns the error induced by evaluating higher order perturbative corrections in an effec- 
tive theory, derived by letting mt — > oo, where mt is the top quark mass. The observation 
that the next-to-leading order (NLO) corrections to the inclusive total cross section are 
approximated at the percent level |6j[7] in this approach has been used as an argument for 
trusting it also at higher orders: next-to-NLO (NNLO) QCD corrections [8-10 have thus 
led to a perturbatively robust prediction for this quantity^ A few years ago, the effective 
theory approach was tested at NNLO by an explicit calculation of the subleading terms in 
1/mt to the total inclusive cross section 20-23 . It was found that these terms have an 
effect of less than 1%. 

However, this result does not allow for a direct generalization to less inclusive quantities. 
Since they depend on several kinematical parameters such as the transverse momentum 
Pt, the rapidity y, or simply phase space cuts, such observables may have a very different 
convergence behavior in 1 /mt than the inclusive cross section. 



Nevertheless, so far the NLO p^- and y-distributions in if+jet-production 24-27 
jet-vetoed Higgs cross section 



the 



28 , as well as the resummation of the logarithmic terms 



for small pi 29-32 are based on the effective theory approach, of course, and so are 



the fully exclusive NNLO partonic Monte Carlo programs for Higgs production in gluon 



fusion 33-35 



Only rather few studies have been aimed at quantitatively testing or going beyond the 
heavy-top approximation in Higgs distributions. Such results are only available at leading 
order (LO) in perturbation theory, for if+rajet with n = 0, 1, 2 [36-40 
one finds that the approximation works rather well for < m t- 



Generally speaking, 



In this paper, we evaluate the 1/mi-effects to -ff+jet production at NLO QCD. We focus 
on the initial states gg and qg which, in the MS scheme, are typically about two orders 
of magnitude larger than all other channels combined. We will show that the radiative 
corrections to the mass effects in the gg channel are remarkably close to those of the heavy- 
top limit. In the sum over all partonic sub-processes, this is deteriorated to some extent 
by the gg-channel which is numerically subleading, however. As a result, we find that 
the kinematical distributions can be calculated with 2-3% accuracy by reweighting the LO 
distributions (including the full top mass dependence) by the differential NLO K-factor 
evaluated in the heavy-top mass limit. 

The remainder of the paper is organized as follows: In Section [2j we briefly specify the 
problem under consideration by introducing the relevant Feynman diagrams; Section [3] 
presents some LO motivation of our study; Section [4] contains the main part of the paper, 
including the results for the pr-integrated cross section as well as for the pi and rapidity 
distributions, all at NLO; our conclusions are given in Section [5] 



1 Effects beyond NNLO and electro- weak corrections have been studied in Refs. 11 ■ 19 , for example 
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(a) 



(b) 



(c) 



Figure 1: Sample Feynman diagrams contributing to the process pp — > ii+jet at 
LO QCD. The graphical notation for the lines is: thick straight = top quark; thin 
straight = light quark q £ {u, d, c, s, b}; spiraled = gluon; dashed = Higgs boson. 



2 Outline of the problem 



In this paper we consider the quantities da/dpT and da/dy in the gluon fusion process, 
where a Higgs boson is produced in association with a jet in hadronic collisions through a 
top-loop mediated gluon-Higgs coupling. Other quark-loop contributions are suppressed 
by their Yukawa coupling and will be neglected. The Higgs' transverse momentum pr 
and its rapidity y are measured relative to the hadronic center-of-mass system. The 
LO contribution to this process is of order a^; it is obtained by convolving the partonic 
subprocesses gg — > Hg, qg — > Hq, qg — )■ H + q, and qq — > Hg (q E {it, d, s, c, 6}), see Fig.[TJ 
with the corresponding parton density functions. At this order of perturbation theory, the 
full dependence on the top quark and Higgs boson mass mn is known, and also parton 



shower effects have been evaluated 39 40) . 

At NLO, the Feynman diagrams can be divided into three groups: the first one is obtained 
by dressing each of the partonic LO processes by a virtual or a real gluon, see Fig.[2](a)-(d), 
for example; the second one by splitting the emitted gluon into a (/(/-pair, see Fig.[2](e). 
The third group is of the form q±q2 — > Hqiq2, where both qi and (/2 run continuously from 
the inital to the final state, and qi, q2 denote quarks or anti-quarks of the first five flavors, 
see Fig.|2](f ) , for example. 



3 Leading order considerations 



Fig. [3] shows the LO result for the cross section 

a( PT >p^)= [ 



A d(J 
dpT -. — 

PT>Pt* PT 



(1) 



as a function of run, divided into the individual partonic sub-processes according to the 
MS-scheme, and keeping the full top mass dependence (solid), the expansion in 1/mt 
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through l/m° (dotted), and through 1/mf (dashed). Unless stated otherwise, we will set 
pS 1 = 30 GeV in this paper; also, we choose the renormalization and factorization scales 
to be /lir = fj,p = mu'i the on-shell top quark mass is set to mt = 172 GeV, and the default 
hadronic center-of-mass energy is y/s = 14 TeV, but we will include exemplary results for 
^/i = 7TeV below. 

The kink in the cross section at mu ~ 2mt in Fig. [3] is due to the top-quark threshold in 
the scattering amplitude. Clearly, this structure cannot be reproduced by an expansion 
in 1 /mt . Note that the shape of the curve is very reminiscent of the total inclusive cross 
section pp — > H + X through gluon fusion which, however, is a 2 — >■ 1 process at LO and 
starts at O(a^); the quantity displayed here is the cross section for H+jet production and 
therefore of 0(o? s ). 

As is obvious from Fig.[3j the l/m^-expansion for the qq channel is significantly worse than 
for the other two channels. This failure of the l/m^-expansion to reproduce the (/(/-channel 



has already been observed for the total inclusive cross section in Refs. 21 , 23 . However, 
the (/(/-channel is about two orders of magnitude smaller than qg which itself is a factor of 
2-3 smaller than gg. Similar observations hold for the other purely quark induced channels 
which enter at higher orders, specifically qq, qq' , and qq' . The conclusion to draw from this 
is that the quark-induced channels constitute a solid, but rather minor limitation of the 
heavy-top limit. Our analysis cannot bring any further insights for these channels, and we 
will disregard them in what follows. One should keep in mind, however, that kinematical 
cuts could enhance the pure quark channels; in this case, results based on the heavy-top 
limit become unreliable. 

The relative deviation of the 1/mt expansion from the exact result (still at LO) is shown 
in Fig.[4](a). The curves are obtained by dividing the integrated cross section as defined in 
Eq. ([I]) when expanded in 1/mt by the same expression when the full top mass dependence 
is kept. The individual plots show this ratio separately for the case when only the gg- 
channel (left) and only the (/(/-channel (center) is taken into account (both in the numerator 
and the denominator of the ratio), and also for the sum of both channels (right). At this 
point, despite our default choice for the center-of-mass energy of ^/s = 14 TeV, we also 
include results for yfs = 7 TeV in order to obtain an impression of the dependence of our 
results on yfs. The corresponding plots are shown in Fig.[4](b). In the following discussion, 
the numbers for yfs = 7 TeV are referred to in brackets; the unbracketed numbers are for 
the default ^fs = 14 TeV. 

For the ^-channel, the heavy-top expansion through 0(l/m|) approximates the exact 
result up to 2% (2.5% for ^/s = 7 TeV) within m H G [100,200] GeV, while the leading 
0(l/mt) term deviates up to 15% (16%) from it. For the (/(/-channel, the approximation 
by the 1/mf -result is not as good: while the l/m° term remains within about 7% (5%) 
of the exact result, the l/m^-term deviates by up to 18% (9.5%). However, since the gg- 
channel is numerically dominant, for the sum of both channels the l/m^-approximation 
agrees with the full result to better than 6% (5%), while the difference between the l/m® 
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and the exact result ranges up to 12% (12.5%). 

We see that the dependence of the 1 / mj-effects on the center-of-mass energy is very weak 
for the dominant ^-channel; for the ^-channel, the l/mj-effects are more pronounced for 
higher ^fs due to the fact that, technically, nit is always assumed to be the largest scale in 



the problem (for a more detailed discussion, see Ref. 20,22,41 ,42]). In order to study the 
quality of the heavy-top approximation, it is therefore sufficient to use our default setting 
= 14TeV. 

The absolute size of the mass effects at LO is used here only as an indicator of how 
far we can expect to be able to trust the heavy-top expansion at NLO. This indicator 
can be considered as a lower limit for the validity range though: In the total inclusive 
cross section, it has proved useful to factor out the LO top mass dependence from the 
perturbative corrections. Our results will show that it is very advantageous to follow this 
strategy also for differential cross sections. 

Turning to more exclusive quantities, Fig. [5] compares the exact result for the Higgs' trans- 
verse momentum distribution to expansions including successively higher orders in l/mt, 
separately for the gg- and the gg-channel. The l/m^-expansion works very well, roughly 
up to pr = nit, as long as one restricts oneself to lower orders in l/m t . At 0{l/m^) and 
beyond, convergence seems to be lost at much lower values of pr- This is due to the region 



of large partonic center-of-mass energy vs: similar to the calculation of Ref. 20-23 , the 
l/mt expansion generates a power behavior in s/m\. At lower orders in l/mt, these terms 
are suppressed by the parton densities. At higher orders, however, they spoil the conver- 
gence behavior of the hadronic cross section. In our NLO analysis, we therefore restrict 
ourselves to comparisons of the l/m\ to the l/m\ terms. 

The behavior of dcr/dpr suggests to try to improve the approximation of the integrated 
cross section a(pT > p™*) by discarding the l/mf-effects above ^ 150 GeV, and taking 
into account only the leading 1 /m® terms for larger values of px • The result is shown 
in Fig.|6j While the effect of this procedure is small in the ^-channel, the contribution 
from > 150 GeV is much more significant in the ^-channel. Even though the numerical 
approximation improves, in particular for the sum gg + qg, this rather large effect indicates 
that the result depends quite strongly on the specific upper cut on pi which has been 
introduced for the l/m\ terms. We conclude that this does not allow for a systematic 
improvement and will not consider it any further. 

Overall, the LO observations are encouraging to study the behavior of the l/mt terms 
at NLO in order to estimate the validity range of the heavy-top limit also for differential 
quantities. 
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Figure 4: Ratio of the integrated LO cross section from Eq. ([I]) when expanded 
through to the exact result, for n = (dotted) and n = 2 (dashed). Left: 

only gg; center: only qg; right: sum of gg and qg; (a) yfs = 14 TeV — (b) 
= 7 TeV. 
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Figure 5: Differential cross section da/dp^ at LO QCD, where pi is the trans- 
verse momentum of the Higgs boson. Solid curve: full mt-dependence in- 
cluded; dotted/dash-dotted/short-dashed/long-dashed: expansion in 1/m™ with 
n = 0/2/4/6. Left, center, and right plot show the gg-, the gg-channel, and their 
sum, respectively. 
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Figure 6: Similar to Fig.JlJ but with an upper pp-cut for the l/m^ -coefficient: 
1/mf -»• l/mj ■ 9(150 -p T/ GeV). 



4 Next-to-leading order results 



4.1 Outline of the calculation 



The most complicated Feynman diagrams are of the two-loop box-type with massless and 
massive (mass mt) internal and one massive external line (mass mjj), see Figs.[2](a) and 
(c), for example. Although not out of reach, the complexity of the corresponding integrals 
is too high for an efficient numerical evaluation. Therefore, to date the NLO corrections 
to this process are only available in an effective theory approach where the top quark is 
integrated out |24| - [27| . The Feynman diagrams then simplify to one-loop level, with an 
effective Higgs-gluon vertex, multiplied by a Wilson coefficient which can be evaluated 



perturbatively 143-47 



The effective theory approach can be seen as the leading term of an expansion for small 
1/mt- The goal of this paper is to go beyond this limit and to study the behavior of 
the next term in this expansion. In Ref. 20 , the relevant one-loop 2 — > 2 and tree-level 



2 — > 3 amplitudes have been obtained through automated asymptotic expansions [48 - 50 
For our purposes, we combine them here with the dipole subtraction terms [51] which we 
have to take into account order by order in 1/mt, of course. The result is a NLO Monte 
Carlo program for ff+parton production in gluon fusion which, in addition to the already 
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available pure heavy-top limit (24] - [26| , also includes the first formally subleading term in 
\/m\. 

We have performed a number of checks on our results. The two most important ones 
are the numerical comparison of the leading terms in 1/mt with the non-resummed part 

The 



of the program HqT 29 30 52] where we find agreement at the sub-percent level, 
amplitudes for the 1/m^-terms have been checked previously by the agreement of the 
inclusive cross section between Ref. 20 and [22]. Their proper implementation into a 



Monte Carlo program is checked by the independence of the numerical results on the 



so-called o-parameter 53,54 which allows to restrict the phase space of the dipole terms. 



4.2 Notation 

We introduce the following notation for the individual terms in our expansions: 

Xe{LO,NLO}, i,je{q,q,g}, (2) 

where I denotes the order of perturbation theory, k the order of the expansion in 1/mt, 
and X the order of the PDFs and the running of a s that have been usecj^J The subscript ij 
denotes the particular partonic channel that was taken into account. If any of the indices 
I, k or ij are absent, it means that these indices are summed over all possible values. For 
example, 

dfT LO = [ d(T (0)]LO ; d(J NLO = [da (0) + d(J (l)]NLO ^ (3) 

are the LO and the NLO differential cross sections with exact mj-dependence, and summed 
over all parton channels (recall, however, that we neglect all qq and qq contributions in 
this paper). 



In order to isolate the individual corrections to the cross section, we define the quantities 



! " da LO {b) 



[K n {b)\ij = 



(4) 



On the right hand side of these definitions, it is understood that da(b) is integrated over all 
kinematical variables except the set b, where we consider b = {pt}, b = {y}, and 6 = (i.e., 
transverse momentum and rapidity distributions, and the integrated cross section with 
Pt > 30GeV). Also, if ij is to be summed over, this applies separately to the numerator 



2 We use the central MSTW2Q08 PDF sets: related uncertainties are not the subject of this paper. See 
Refs. [551156] , however. 
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and the denominator in [K n (b)]ij. The ratio R k \b) allows for a direct comparison of the 
perturbative (index I) and the mass effects (index k). The quantity K(b), on the other 
hand, shows the influence of the mass terms on the perturbative correction factor. For 
example, Kq is the NLO K-factor in the heavy-top limit which - in the case of the total 
inclusive cross section - has been found to approximate the exact NLO K-factor extremely 
well. Using the l/mj expansion, we will study whether this observation can be expected 
to carry over also to differential quantities. 



4.3 Inclusive Higgs plus jet production 

The first observable we study is the integrated cross section for Higgs+jet production, 
defined in Eq. @. Fig.[7] compares the NLO perturbative corrections to the mass effects 
at LO and at NLO, split into the two numerically dominant sub-channels gg (dotted) and 
qg (dashed) , as well as for the sum of both channels (solid) . The size ofiif (=4 O) (0)) 
is mostly determined by the reduced value of a s when going from LO to NLO parton 
densities. R*q\ on the other hand, reflects the well-known largeness of the perturbative 
effects to the gluon fusion cross section. Considering the fact that JSp includes a factor 
a s /ir relative to R , it is remarkable that they are both almost equally large. Note also 
that both and R^ depend only very weakly on the Higgs mass rriH- 

The same feature holds for the mass effects, shown in the lower two plots, separately 
for the LO (left) and the NLO (right) coefficients. Note also that there is a cancellation 
between the gg and the qg channels, although much less pronounced at NLO than at LO. 
The overall mass effects between m# = 100 GeV and 200 GeV range from —2% to 6% for 
the LO, and from 2% to 8% for the NLO coefficient. They are thus much smaller than the 
perturbative effects. As expected, the mass effects decrease for smaller Higgs masses. 

Concerning the K-factor, for the total inclusive cross section it has been found to depend 
only very weakly on the top quark mass [6j[20|[22]. The product of the K-factor with the 



exact LO cross section is thus an excellent approximation of the higher order cross section. 

Let us study the extent to which we can draw a similar conclusion for the cross section 
with a lower cut. In Fig.|8j we compare the result for the K-factor including mass terms, 
K2 (= K2{%), cf. Eq. Q) to the pure heavy-top limit. Again, we consider separately the 
channels gg and qg, as well as their sum. Note however, that according to the definition 
in Eq. Q, [-fT n ]jj really only refers to the ij channel, both in the numerator and the 
denominator. Therefore, K n 7^ [K n ]gg + [K n ] qg . 

The agreement between Kq and K2 for the gg channel is truly remarkable; for the qg 
channel, we find 5-10% difference, but due to the numerical dominance of gg, the overall 
agreement between Kq and K2 is around 3%. 
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Figure 7: Relative size of the perturbative and the mass effects on the integrated 
cross section, = i?j^(0), see Eq. |l]) and Eq. (4). Upper row: NLO effects 
arising solely from the PDFs and the running of a s (left), and from the NLO 
perturbative coefficient in the cross section (right). Lower row: mass effects from 
the LO (left) and the NLO (right) perturbative coefficient. Dotted: only qg; 
dashed: only gg; solid: sum of gg and qg. 
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Figure 8: K-factors as denned in Eq. Q, for the integrated cross section, i.e. 
K n = K n ($). Left/center/right plot: [K n }gg/[K n } qg /[K n ] gg+qg . Dotted/dashed: 
n = 0/2. The dots show the results of our calculation; the lines have been intro- 
duced to guide the eye. The deviation between the dots and the lines indicates 
our numerical error. 
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4.4 Transverse momentum distribution 



Fig. [9] shows the pr-dependent ratio R(pr) (= R{{pt}), cf. Eq. Q), in analogy to Fig. [7} 
The qualitative features of the individual corrections for this differential quantity are very 
similar to the integrated ones. An observation that deserves to be pointed out is the 
stunning similarity of the plots for the mass effects at LO and NLO. 

Both at LO and NLO, the mass terms in the gg-channel are very small (~ 2%) and almost 
independent of pt, even up to pr = 300 GeV. The qg-channel behaves worse, but its mass 
terms still do not amount to more than 6% below pi = 150 GeV; at LO, they reach up 
to almost 50% at = 300 GeV though. The NLO mass terms are only slightly smaller. 
The sum of gg and qg, however, remains below 35% at NLO for pr < 300 GeV; below 
Pt = 150 GeV, they amount to not more than 3%. The mass effects reduce the absolute 
value of the cross section in the qg-channel significantly for > 150 GeV, which is also 
visible at LO in Fig. [5} while it just slightly affects the sum of both channels. 

The K-factors 
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The px-dependent K-factors K n (pr) (= K n ({p^)) are shown in Fig. 
including leading and subleading mass terms, Kq and K2, are almost identical in the gg- 
channel. For the gg-channel, on the other hand, the QCD corrections to the subleading 
mass terms behave very differently to the leading terms in 1/mt once px > 150 GeV. In 
the sum of both channels, the difference remains below 3% for p^ < 150 GeV, and reaches 
10% at p T = 300 GeV. 

In conclusion, the behaviour of K2 with respect to Kq suggests that, also for the pt- 
distribution, the QCD corrections can be safely calculated in the heavy-top limit; the 
accuracy remains within 2% (10%) below p^ = 150 GeV (pr = 300 GeV). The absolute 
distributions, however, should be calculated at LO using the full top-mass dependence, 
and then reweighted by these QCD corrections. 



4.5 Rapidity distribution 

Fig. [ll| shows the y-dependent ratios R(y) = R({y}). The perturbative effects, -Rg ^ and 
Rjp , are qualitatively very similar to the quantities discussed before. The mass effects 
are generally very small over the full y-range. At LO, there is a significant cancellation 
between the gg- and the qg- channel, so that the sum remains below 1% almost everywhere, 
even though the individual channels reach up to 4%. At NLO, the mass effects in the qg- 
channel are very small, but due to the numerical dominance of the purely gluon induced 
contributions, the overall effect reaches up to 4%. 



The K-factors, shown in Fig. 12, display a similar behavior as in the previous observables: 
the corrections in the gg-channel are practically the same in the l/m^- and the 1/mf- 
terms, the gg-channel shows some difference, but in the sum of both channels, the K2 
is approximated by Kq to within about 3%. Apparently, the bad convergence of the qg- 
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Figure 9: Similar to Fig.JTJ but for the differential cross section da/dpT, here, 
4° = 4°(-br})- 



channel for px > 150 GeV as observed in section 4A affects pr-integrated quantities only 
at the percent level, see also Fig.|8l 



5 Conclusions 



In this paper, the quality of the heavy-top limit in the gluon fusion cross section has been 
studied. Subleading terms in l/m t have been calculated for the .ff+jet cross section at 
NLO QCD, and their effects on the Higgs' transverse momentum and rapidity distribution 
have been evaluated. 

We found that, similar to the leading terms in l/mt, the perturbative corrections on the 
subleading terms are of order one. In fact, the perturbative effects on the mass corrections 
in the gg-channel are remarkably similar to those on the leading mass terms, both for 
the pt- and the y-distribution, as well as for the integrated cross section of Eq. Q. The 
NLO K-factors with and without mass terms are therefore almost identical for this channel 
alone. Including the gg-channel spoils this similarity to some extent, but we still claim 
that the procedure of correcting the full LO prediction (including top mass effect) by the 
K-factor as evaluated in the heavy-top limit provides an excellent approximation to the 
full NLO result, valid at the 2-3% level for pi < 150 GeV and for px-integrated quantities. 
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Figure 10: Similar to Fig.[8j but for the differential cross section dcr/dpr; here, 
K n = K n ({p T }). 



We have checked that this result holds for Higgs masses below 2mt- The accuracy is thus 
better than the current uncertainty on the cross section due to its dependence on the PDFs 
and due to missing higher order QCD corrections. 
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